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An alternative equation, resumming of the In 1/x terms for the polar- 

^—v , ized nonsinglet structure function g± at small x is presented. Construction 

~~ — of the GLAP-like formula for the auxiliary function, corresponding to the 

i-£h J g^ s a t rescaled Q 2 variable is shown. Predictions of this approach for the 

*~" | i g^ s function at small x in a case of a flat as well as a dynamical input are 

f*) |' given. The role of the fixed coupling constant and the running one is also 

QJ , discussed. 
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» . 1. Introduction 

Our knowledge about the structure functions of the nucleon is still in- 
complete because of lack of understanding how these structure functions 
behave in the small Bjorken x region. Neither present experimental data 
nor the theoretical QCD description give a full and unique picture of an 
exact shape of the quark and gluon distributions in the nucleon. Perturba- 
tive QCD analysis, based on the GLAP evolution equations [1] is in a good 
agreement with experimental measurements. This agreement concerns un- 
polarized [2], [3] and polarized [4] structure functions of the nucleon within 
NLO approximation in the large and moderatory small Bjorken x region. 
Unfortunately, practically lack of experimental data in the small x region 
(x < 1CP 3 ) makes the satisfactory determination of the Bjorken and Ellis- 
Jaffe sum rules [5] impossible. This causes e.g. that the question "how is 
the spin of the nucleon made of partons?" remains still open. From recent 
papers [6], [7] we know, that the small x behaviour of both unpolarized and 
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polarized structure functions is controlled by the double logarithmic terms 
a s In 1/x. However, this singular behaviour of the structure functions at 
low x is better visible in the spin dependent case. For the unpolarized, non- 
singlet structure function Ff s (x, Q 2 ) = F%(x,Q 2 ) - F?(x,Q 2 ) the QCD 
singular behaviour at small x is overridden by the leading Regge contribu- 
tion [8]. Next in the unpolarized, singlet case, the structure function at low 
x is driven by BFKL pomeron [9] because gluons play the dominant role. 
Thus the growth of structure functions of the nucleon, governed by leading 
double logarithmic terms a" In n (x) becomes best visible for spin dependent 
functions. Therefore the polarized structure functions of the nucleon may 
be a sensitive test of the perturbative QCD analyses in the small x region. 
The double logarithmic In x effects go beyond the standard LO or even 
NLO QCD evolution of parton distributions and correspond to the ladder 
diagrams with quark and gluon exchanges along the ladder. One has also to 
take into account nonladder diagrams but in the nonsinglet case they may 
be neglected as nonleading [6], [7]. Thus the nonsinglet, polarized struc- 
ture function gi S (x,Q 2 ) = gi(x,Q 2 ) — g"i(x,Q 2 ) is a convenient function 
both for QCD analyses (because of its simplicity) and future experimental 
tests at HERA [10] concerning the determination of the Bjorken sum rule. 
Theoretical predictions for gi at small x, incorporating the double In x 
effects have been presented in [7], [11], [12], [13]. In these papers the per- 
turbative QCD analysis is based on the unintegrated spin dependent quark 
distributions f(x,k 2 ). It means that the sum of double logarithmic In x 
terms is represented by a appropriate integral equation for the unintegrated 
structure function f(x,k 2 ). 

In this paper we present an alternative approach of the double In x 
resummation for g^ at small x. Our formalism is based on the usual 
quark distribution functions (and not on f(x, k 2 ) function) and moreover 
generates the In x terms via GLAP-like evolution equation for the rescaled 
quark transverse momentum squared fi 2 = k 2 jx. The purpose of this paper 
is to compare these two methods of the double In x resummation for the 
polarized nonsinglet structure function g± s at low x. In the next section 
we briefly recall the origin of the double logarithmic In x effects at low 
x, incorporating the evolution equation based on the unintegrated function 
f(x,k 2 ). In section 3 we introduce alternative formalism, in which after 
rescaling the kinematic variable fj, 2 = k 2 jx we get GLAP-like evolution 
equation in fi 2 . This equation for the polarized quark distributions (and 
hence for the g^ function) generates the double logarithmic In x terms. 
Section 4 contains numerical results for the spin dependent nonsinglet struc- 
ture function g^ (x, Q 2 ) in our approach. We compare the both mentioned 
above methods and also compare their predictions for g^ with SMC 1997 
small x data. Finally in section 5 we summarize our results. 



2. Double logarithmic In x resummation for the nonsinglet, 
polarized structure function g™ s (x,Q 2 ) using the 
unintegrated function f(x,k 2 ). 

It has been lately noticed [6], [7] that the spin dependent structure func- 
tion gi in the small x region is dominated by In (1/x) terms. This singular 
behaviour, implied by QCD is for the polarized structure functions the lead- 
ing one. Comparatively, for unpolarized, nonsinglet structure functions of 
the nucleon, the QCD evolution behaviour at small x is screened by the 
leading Regge contribution. The Regge theory [14], which concerns the 
Regge limit: x — ► predicts the following behaviour of parton distributions 
at small x and Q 2 < 1 GeV 2 : 

xS ~ const (Pomeron) 
Qns ~ x' - 5 (Reggeon A 2 : p - u); 
AS, Aq NS ~ x°-hx - 5 {Reggeon A x ) (2.1) 

where S, qNS-, AS, Agjvs denote respectively singlet unpolarized, nonsinglet 
unpolarized, singlet polarized, nonsinglet polarized quark distributions. The 
shape of all spin dependent distributions AS, Aq^s is mostly governed by 
QCD evolutions with dominating In x terms at small x. These In x contri- 
butions correspond to the ladder diagram with quark and gluon exchanges 
along the ladder - c/Fig.l. In contrast to the singlet polarized function, 
for the nonsinglet one the contribution of nonladder diagrams is negligible. 
Thus examing the polarized, nonsinglet structure function g^ s (x,Q 2 ), we 
should consider only mentioned above ladder diagrams. The nonsinglet part 
of the spin dependent structure function has a form: 

g? s (x,Q 2 )=g p 1 (x,Q 2 )-g?(x,Q 2 ) (2.2) 

where g± and gf are spin dependent structure functions of proton and neu- 
tron respectively. Let us remind the meaning of g\. In the Bjorken limit 

1 

=u,d,s,.. 

Aft(x) = qi+{x) - qi-(x) (2.4) 

where e$ is a charge of the i- flavour quark, </j+(x) (qi-(x)) is the density 

distribution function of the i-quark with the spin parallel (antiparallel) to 

the parent nucleon. Function gi(x,Q 2 ) is connected with the helicity of 

the nucleon (i.e. spin projection on the momentum direction). Thus the 

integral 

l 

(Aft) = J A qi {x)dx (2.5) 

o 



9i(x) = - E e ? A %0*0 ( 2 - 3 ) 
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Fig. 1. A ladder diagram generating double logarithmic In (1/x) terms in the non- 
singlet spin structure function g\. 



is simply a part of the nucleoli helicity, carried by a quark of i-flavour 
(i=u,d,s,..). Polarized distribution functions of quarks are denned as: 



Aq = Aq va i + Aq s 



(2.6) 



Finally: 



<?f 5 = -{Au val 
b 



Ad 



val j 



-(Aw -Ad) 





(2.7) 



Solutions of the equation for the unintegrated polarized nonsinglet structure 
function / (x, Q 2 ), which gives the In x resummation are presented in [7], 
[11], [12], [13]. In x as the only contribution at small x behaviour of g^ 
is examined in [7] while in [11], [12], [13] the unified description of g^ s 
incorporating both Gribov-Lipatov-Altarelli-Parisi (GLAP) evolution and 
In 2 x effects is presented. These all approaches are based on the unintegrated 
distribution function f(x,Q 2 ), which is related to the gi(x,Q 2 ) via 



Q 2 (l/x-l) 



gi(x,Q 2 



9? ) (x) + 



I.-2 



dk 2 , ( , fc\ 



where 



I.-2 



dk 2 



g?\x)=j^f(x,k 2 ) (2.9) 

o 

Let us recall from [11] that the resummation of the double logarithmic 
terms In x in the limit of a very small x (x — ► 0) is given in the case of the 
polarized, nonsinglet structure function by: 

1 k 2 /z 

fNs(x,k 2 ) = fSl(x,k 2 ) + ^lj^ J ^AP(J(,)fe(^' 2 ) 

X 1.2 

(2.10) 
The source of the double logarithmic terms In x in g\(x, Q 2 ) is the double 
integration in the formula for function f(x,k 2 ): 

1 k 2 /z 

,, , 9n a s f dz f dk' 2 

'<*•*> ~s/t / -w (2 - u) 

X 1.2 

K Q 

where the upper limit in the integral over the transverse momentum k' 2 
is z-dependent (= k 2 /z). Thus, double logarithmic terms come from the 
integration over the longitudinal momentum fraction z together with the 
integration over k' 2 with ^-dependent upper limit: 

f(x,k 2 ) ~ In 2 (l/x)= In 2 x (2.12) 

The analytical solution of (2.10) in the case of the fixed coupling constant 
a s [12] shows the singular small x behaviour of the polarized, nonsinglet 
structure function g\ i.e.: 

g? s (x,Q 2 )~x~ x ^; (2.13) 



\ns = 2 ] /^AP£>(x); (2.14) 

where APq q (x) is the splitting function and in the limit x — ► is equal 
to 4/3. Hence for the fixed a s =0.18, Xns = 0.39 and as it has been al- 
ready mentioned above, the singular small x shape of g^ , implied by QCD 
dominates the REGGE behaviour: 

REGGE : g? s (x,Q 2 ) ~ x ~ aNS ^ a NS (0) < (2.15) 

QCD: g? s (x,Q 2 ) ~ x~ Ajvs X NS ~ 0.5 (2.16) 



In the case of the running coupling constant, numerical results for g^ , 
based on the Regge flat input parametrization gi (x,Qq = 1 GeV 2 ) are 
given in [7], [11], [12]. In [13] the same method is used for dynamical 
input parametrization, where in contrast to Regge flat parametrizations, 
gi (x, Qq = 1 GeV ) is singular for small x. The generation of the double 
logarithmic In x terms is also possible via GLAP-like evolution equation 
for gi function with the rescaled transverse momentum squared. This al- 
ternative method is presented in the next section. 



3. GLAP-like evolution equation for g± , generating In x terms 

at low x. 

We study the origin of the double logarithmic In x terms in g^ s {x, Q 2 ) 
using a simple rescaling of Q 2 variable: Q 2 — ► li 2 = Q 2 /x, what leads to the 
GLAP-like equation for g± s with evolution in a new scale li? . We focus on 
the small x region so our initial equation for further investigations is that 
which contains only dominating In x part. First let us consider the case 
with a fixed coupling d s = const, where 

a s = — (3.1) 

Thus the starting equation is: 

1 k 2 /z 

fNs(x,k 2 ) = f^l(x,k 2 )+a s J^ J ^f NS ^,k' 2 ) (3.2) 



X u2 

K Q 



After a simple substitution: 



h 2 t k' 2 

fj, = — ; x = -; (i =— (3.3) 



X Z X 

the equation (3.2) takes a form: 



f NS (x,xn 2 ) = f$' s (x,xfj?) + a s / — / —pffNsi-,-^ 2 ) (3.4) 

J Z J LL Z Z 



k 0x 



After introduction an auxiliary function </?(x,/x 2 ): 

tp(x, /i 2 ) = /ns(x, xix 2 ) (3.5) 



and applying Heaviside's G function: 



e«)Hj tt<l < 3 ' 6 > 



we get: 



i . m 2 



rts,M 2 ) = ¥> (0 W)+*- I d - l^e^-^-M-.M") (3-7) 

J Z J /i' z X z 

a; 1.2 

or 

1 , ^ , ,i 
i 2n (o), 2 ,, - [ dz [ d/i" x , 2 
¥>e(a, n ) = <f @ {x,fi )+a s I — —^-(f e (-,fi ) (3.8) 

Z J jji z 

X 

where 



IE 1.2 



^e(x,/z 2 ) = e( M 2 -^Mx,M 2 ) (3.9) 

x 

The equation (3.8) has an exact form of the GLAP Q 2 evolution formula for 
the unintegrated distribution function f(x,Q 2 ). The mentioned GLAP evo- 
lution equation for polarized nonsinglet quark distributions Ap (and hence 
for g^ s function too) in the small x region is given by: 

l 
dAp(x,Q) _ f dz x 



dlnQ 2 



a s [ —Ap{-,Q 2 ) (3.10) 

J z z 



Relation between Ap(x, Q 2 ) and the unintegrated distribution f p (x,Q 2 ) is 
as usual: 

AC*) -^ (3-"> 

what implies: 

O 2 

r do' 2 
Ap(x, Q 2 ) = A Po (x) + J ^rf P (x, Q' 2 ) (3-12) 

Ql 
where Apq{x) is a nonperturbative part of Ap: 



Q f°dO' 2 
A P0 (x) = j -^f p (x,Q' 2 ) (3.13) 



Q 





and Qq = 1 GeV is the low scale of perturbative QCD. Hence the evolution 
equation (3.10) written for the unintegrated distribution function f p (x,Q 2 ) 
takes a form: 

f P (x,Q 2 ) = ff\x,Q 2 ) + a s J ^ J B^f^Qp) (3.14) 

x Z Ql 



and 



l 

dz . ,x. 



4°)(x,Q 2 ) = a s /-Apo(-) (3.15) 

J z z 

x 

One can see from (3.8) and (3.14) that the double In x resummation equa- 
tion written for (p®(x,fj?) function is a GLAP fj, 2 evolution equation. The 
auxiliary function (Pq(x,/j, 2 ) may be, similarly as in (3.11), represented by 
an integrated function n(x,/x 2 ): 



2 du(x,fi 2 ) 
Mx,/i) = -gj^- (3-16) 



and conversely: 



f 



V 2 ,. , . .,2, 



u(x,n ) = u (x) + / —^-(f Q {x,n ) (3.17) 

Thus the equation (3.8), generating double In 2 x effects can be rewritten as: 

l 
du(x,fi 2 ) _ f dz x 2 

- a s / —u{-,n ) (3.18) 



cHn^ 2 J z z 



X 



Relation between the nonsinglet polarized structure function g±(x,Q 2 ) 
and the auxiliary function u(x,fj?) is as follows: 

gi S (x, Q 2 = xfi 2 ) = u(x, f?) (3.19) 

In this way the problem of producing the In x terms for g^ s (x, Q 2 ) in the 
small x region via equation (2.10) has been reduced to the GLAP evolution 
to the momentum scale /i 2 = Q 2 /x. It is not astonishing: appearing of 
the new evolution scale Q 2 /x has its origin in the upper limit k 2 / z of the 
integration over the transverse momentum in (2.10). This logarithmic inte- 
gration over the transverse momentum up to the z— dependent limit k 2 /z 



together with the logarithmic integration over the longitudinal momentum 
fraction z give double logarithmic ln 2 x terms. The mechanism of appearing 
of the ln 2 x effects in g^ s from GLAP-like equations (3. 18)- (3. 19) is well 
visible just in a case of the fixed coupling constant a s = const. Then the 
eq.(3.18) can be solved analytically. Using standard Mellin's method one 
can get the solution of eq.(3.18) in the form (see Appendix A): 



,2 . 



^ (a s ln±Znf7) fc 



«0^ 2 )~£ — mr^ (3-20) 

k=0 



k\k\ 



and hence: 



^(.,«ii^4)~£ ' • ; kl ¥ (3.2D 

X k=Q K.K. 

what gives approximately the leading small x behaviour: 

<7f 5 (x,Q 2 )~af- 2 v^ (3.22) 

Taking into account parton interactions through the introduction of the 
running coupling constant one can get from (3.18) an equation which incor- 
porates the running couplings effects a s — ► a s (Q 2 ): 

du(x,n 2 ) _. 2 \ f dz ,x 



<91n/i 2 



a s ( Xf i 2 ) f- u (-^ 2 ) (3.23) 

J z z 



However more justified theoretically seems the introduction of the running 
coupling by the substitution a s — ► a s {Q 2 / z), what gives: 

°<*,H 2 )_ ^(^)u(-V) (3.24) 



d In /i 2 J z z z' 



X 



Our numerical analysis presented in the next section contains the both above 
"running coupling" prescriptions and the constant a s case as well. 

4. Numerical predictions for g™ s (x,Q 2 ) based on the GLAP-like 
equation, resumming the In x terms. 

We solve numerically equations (3.18), (3.23) and (3.24) which are GLAP 
/x 2 evolution equations for the auxiliary function u(x,/i 2 ). The relation 
between u(x,/x 2 ) and the physical polarized nonsinglet structure function 
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gf (x,Q 2 ) (2.7) is given by (3.19). In this way one can get the small x 
behaviour of gi S governed by the double logarithmic In x effects. Our 
predictions we compare with those, received in the unintegrated f(x,Q 2 ) 
approach (2.10) and described in [7]. Solving the GLAP equations (3.18), 
(3.23) and (3.24) one should have an input parametrization of the u function 
at the low scale fc 2 . Because u{x,Q 2 ) has a meaning of the physical g± 
function for the rescaled Q 2 variable as it is shown in (3.19), one can also 
write: 

u(x,k 2 )=g? s (x,xk 2 ) (4.1) 

The low scale kg introduced in (3.7) is equal to the usually used QCD cut-off 
parameter /c 2 = 1 GeV 2 . The really dependence on z/x of the lower limit 
in the integration (3.4): 

Q 2 = Z -k 2 (4.2) 

became "shifted" to the definition of ipo(x,n 2 ) via (3.9). It means that 
for running coupling cases one should take into account in the evolution 
equations (3.23) and (3.24) the cut-off factor 0(/i 2 — k^/x). Otherwise, the 
running variable x/j? in the coupling a s becomes less than the low scale of 
QCD evolution k 2 = 1 GeV what is definitely incorrect. There is no such 
a constraint in the fixed coupling constant. Taking into account that for 
perturbative QCD analysis the low cut-off parameter k 2 , = 1 GeV we put 
this value for the input scale of the u function too. This implies that for 
small x the input parametrizations u(x, kg) corresponds to the gi S at the 
very low scale xkg. We assume, that below the kg = 1 GeV 2 the behaviour 
of the quark distributions is the same as at fc 2 . Therefore we apply the 
standard parametrizations of the valence quarks (and hence of gf ) for the 
auxiliary function u: 

u(x, k 2 ) = g? s (x, k 2 ); k 2 = 1 GeV 2 (4.3) 

There are two basic kinds of input parametrizations of gi S (x, k 2 ,): the Regge 
one, which is flat at small x and the singular one, which behaves like x~ a (a ~ 
0.3) at small x. In our numerical calculations we use two different input 
parametrizations: the Regge one, which is given by 

REGGE : u(x, jfcg) = -g A (l - xf = 0.838(1 - x) 3 (4.4) 

where g A = 1-257 is the axial vector coupling and the dynamical input 
GRSV [15]: 

GRSV : u(x, kl) = 0.327x-°- 267 

x (1 - 0.583x - 175 + 1.723x + 3.436a; 3 / 2 ) (1 - x) 3 ' 486 

\4.215 



+ 0.027x"°- 624 (l + 1.195x - 529 + 6.164x + 2.726x 3/2 )(l - .r) 



(4.5) 
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For details about these parametrizations see also [13]. In all calculations 
h-QCD = 232 MeV. Our numerical results are presented in Figs. 2-5. In 
Fig. 2 the predictions for g± s at small x, based on the eq. (3.18) for Regge 
and GRSV inputs respectively are shown. We use two different value of 
fixed coupling: a s = 0.18 and a s = 0.12. The input parametrizations are 
also plotted. In Fig. 3 we confront the fixed coupling results of eq.(3.18) 
for a s = 0.18 with those, based on eqs. (3.23)-(3.24), taking into account 
running coupling effects. Two groups of lines correspond to different inputs: 
Regge and GRSV. The effective slope Xjys determined from (2.16) as: 



Ajvs 



dlng{ 
3 In ■ 



NS 



(4.6) 



for both input parametrizations and all mentioned above a s cases: a s = 
const = 0.18, a s (Q 2 ), a s (Q 2 /z) is presented in Fig.4. The comparison 
of small x predictions for g^ s based on the GLAP-like In x approach with 
those based on the unintegrated function f(x, Q 2 ) and the eq.(2.10) is shown 
in Fig. 5. Both parametrizations are used. The running coupling effects 
a s (Q 2 ) are included. We also plot few small x recent SMC 1997 data [4]. 
In all plots Q 2 = 10 GeV 2 . 




Fig. 2. The nonsinglct spin structure function of the proton g^ (x, Q 2 ) received 
from (3.18). Solid lines correspond to Regge input, dashed lines - to GRSV input. 
In each group of lines (Regge or GRSV) the lowest (for small x) line is the input at 

the middle one concerns the fixed coupling a s = 0.12 and the upper 

= 0.18. Evolution scale Q 2 = 10 GeV 2 . 



fc 2 = 1 GeV 
one is for a s 
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Fig. 3. The nonsinglet spin structure function of the proton g^ s (x, Q 2 = 10 GeV 2 ). 
Solid lines correspond to Rcggc input, dashed lines - to GRSV input. In each group 
of plots (Regge or GRSV) the lowest line concerns the running coupling a s (Q 2 /z) 
case (3.24), the middle one corresponds to a s (Q 2 ) case (3.23) while the upper plot 
is for fixed coupling a s = const = 0.18 (3.18). 




Fig. 4. Slope Ajvs defined in (4.6) for different inputs (Regge - solid, GRSV - 
dashed) and different a s . For each group of lines (Regge or GRSV) the upper plot 
corresponds to fixed a s = 0.18, in the middle lies line for the running a s (Q 2 ), and 
the lowest plot is for the "more running" case: a s (Q 2 / ' z). Q 2 = 10 GeV 2 . 
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Fig. 5. Comparison of g± s predictions based on GLAP-like eq. (3.23) - solid with 
that based on fix, Q 2 ) approach (2.10) - dashed at Q 2 = 10 GeV 2 . Triangles show 
the recent small x SMC data 1997 [4]. 



From Figs. 2- 5 one can read the following conclusions: 

l.The In x resummation gives steep growth of the structure function in 
the small x region. It is well known effect [6], [7] that for x < 10 -2 the In x 
terms dominate over the LO (or even NLO) evolution. 

2. The growth of the structure function in the small x region is of course 
much steeper for the dynamical parametrization than for the flat one. Sin- 
gular inputs at fcg intensify the QCD evolution effects while in the case of 
flat parametrizations the singular small x behaviour of structure functions 
is completely and only generated by QCD evolution. 

3. The effective slope X]\rs (4-6) is contained between 0.14 for the non- 
singular input and the running coupling a s (Q 2 / z) and 0.6 for the singular 
input and the fixed coupling a s = 0.18. This shows apart from the above 
conclusion in the point 2 that for the running coupling constant a s (Q 2 ), Xns 
is smaller than for the fixed one. Besides, introduction of the a s = a s (Q 2 j z) 
causes the slope Xns still smaller: 

XNs(a s = const = 0.18) > X NS (a s = a s (Q 2 )) > X NS {a s = a s (Q 2 /z)) 

(4.7) 
For comparison, the Regge theory predicts in the small x region ola x (0) < 
and <7i behaves as x~ aA ^ ' ~ const. 

4. Comparing our results with those, based on the unintegrated function 
f(x,Q 2 ) and the equation (2.10) one can see agreement of these two ap- 
proaches. It is not astonishing because the both methods concern in fact 
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the same problem: resummation of the double logarithmic terms ln 2 x for 
the g^ s function in the small x region. Our approach based on the GLAP- 
like equation (3.18) enables to keep the Constance upper scale of evolution 
/i 2 instead of Q 2 /z like in (2.10). In this way this approach become similar 
to the standard GLAP Q 2 evolution and we can adopt the methods used in 
it to get new - In x effects. 

5. SMC recent data [4] for very small x (3 points shown in Fig. 5) confirm 
the growth of the polarized nonsinglet structure function g^ s in the small 
x region. It seems that the singular inputs enable the better agreement of 
theoretical predictions with experimental data. 



5. Summary and conclusions 

In this paper we have presented the double logarithmic In x terms re- 
summation for the spin dependent nonsinglet structure function g± . In our 
approach we have solved the GLAP-like equation for the auxiliary function 
u{x,n 2 ) which corresponds to the physical function g^ s (x,Q 2 = x/x 2 ) at 
the rescaled jx 2 variable. Our calculations have been performed for the sim- 
ple nonsingular Regge parametrization and for the dynamical one as well. 
Apart from the fixed coupling constant case, the running coupling effects 
have been taken into account. Besides the effective slope, controlling the 
small x behaviour of g± has been estimated. Its value has been found lies 
between 0.14 for the nonsingular input and the running coupling a s (Q 2 j z) 
and 0.6 for the singular input and the fixed coupling a s = 0.18. We found 
that the In x effects govern the small x increase of the structure function 
gi S (x, Q 2 ). This growth is larger in the case of the fixed coupling constant 
than for the running one and of course for the singular input g^ s (x, k^) than 
for the flat (e.g.Regge) one. The equation we have considered is not appli- 
cable for the large x (x > 10 -2 ) region. Our formalism is however correct 
in the very interesting small x region. Presented results confirm that the 
In x effects are very significant for x < 10~ 2 . The spin dependent structure 
functions of the nucleon are a sensitive test of the perturbative QCD anal- 
yses in the small x region. However practically lack of experimental data in 
the very low x region (x < 10~ 3 ) causes the satisfactory verification of the 
theoretical QCD predictions in this region impossible. Also the predictions 
incorporating the double logarithmic In x effects in g^ are still awaiting 
for their crucial probe. 
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Appendix A 

Analytical solution of the GLAP-like evolution equation for g^ s generating 
double logarithmic In x effects at small x. 

We solve the equation (3.18) with the fixed coupling d s = 0.038 using 
the standard Mellin method. The Mellin transformation defines the n— 
moment of u(x,/i 2 ) function by: 

i 
u n (fi 2 )= j x n - l u(x,ii 2 )dx (A.l) 



Thus in moment space the evolution equation for small x (3.18) is simply 
given by: 

d«»V) _ a. 2 



dlnfj, 2 n 
The solution of (A. 2) is straightforward: 



"V)=<[p] : (A.3) 



where Uq is the n— moment of the input u(x, Uq) function: 



i 

ra-l„./„. ;„2n 



u%= / x n - L u(x, jfcg)dx (A.4) 



For the nonsingular input parametrization 

u(x, fe ) ~ const (A. 5) 

Uq has a form: 

n const 
«S = — — (A.6) 

Employing the inverse Mellin transformation: 

C+IOD 

u(x,n 2 ) = — I x~ n u n (\x 2 )dn (A.7) 

2iri J 

c—ioo 

which gives 

c+ioo 

2 const f x~ n rn 2 i^ 
u[x,n ) - — J -j dn (A. 8) 



2ni J n 

c—ioo 
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One can get the solution u(x,fx 2 ) of the form: 



u(x, u 2 ) = const J2 ,.,, V (A.9) 



<x> (a K ln^ln 1 U) k 



A-=0 klkl 



Going back to the function g± (x, Q 2 = x/x 2 ), one can find from (A.9) 
through (3.19) the following expression: 



O 2 / 1 O 2 

5 f s (x,Q 2 ) = u(x,— ) = const J (2Ja s ln-ln-^- 1 ) (A.10) 






where Jo(y) denotes modified Bessel function: 

oo (y\2k 

This In x terms resummation (A. 10) gives in the small x region the effective 
behaviour of gi S : 

gf s (x,Q 2 ) = x- XNS (A.12) 

where 

X NS ~ 2sfbT s (A.13) 

For the fixed a s = 0.038 the effective slope Xjys is equal to 0.39. 
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